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Abstract 

We prove that any potential symmetry of a system of evolution 
equations reduces to a Lie symmetry through a nonlocal transforma- 
tion of variables. Based on this fact is our method of group classifica- 
tion of potential symmetries of systems of evolution equations having 
non-trivial Lie symmetry. Next, we modify the above method to gener- 
ate more general nonlocal symmetries, which yields a purely algebraic 
approach to classifying nonlocal symmetries of evolution type systems. 
Several examples are considered. 

1 Introduction 

Lie symmetries and their various generalizations have become an insepara- 
ble part of the modern physical description of wide range of phenomena of 
nature from quantum physics to hydrodynamics. Such success of a purely 
mathematical theory of continuous groups developed by Sophus Lie in 19th 
century p!] is explained by the remarkable fact that the overwhelming ma- 
jority of mathematical models of physical, chemical and biological processes 
possess nontrivial Lie symmetry. 

One can even argue that this very property, invariance under Lie symme- 
tries, distinguish the popular models of mathematical and theoretical physics 
from a continuum of possible models in the form of differential or integral 
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equations (see, e.g., [21 Based on this observation is the symmetry se- 
lection principle stating that if an equation describing some physical process 
contains arbitrary elements, then the latter should be so chosen that the 
resulting model possesses the highest possible symmetry. In this sense Lie 
theory effectively predicts which equation is the best candidate to serve as a 
mathematical model of a specific physical, chemical or biological process. 

The process of choosing from a prescribed class of differential equations 
those enjoying the highest Lie symmetry is called group classification. In the 
case when non-Lie symmetries are involved, the more general term, symmetry 
classification, is used. 

In this paper we study symmetries of systems of evolution equations in 
one spatial variable 



where u = {u^(t, x),U2(t, x), . . . , u"^(t, x)}, Uj+i = dui/dx, n > 2, m > 2. 

There is a lot papers devoted to group classification of different subclasses 
of the class of partial differential equations (PDEs) of the form ([1]) (see, e.g. 
[l]-[7] and the references therein). The major tool utilized in these studies is 
the infinitesimal Lie approach enabling to reduce the problem of exhaustive 
description of Lie transformation groups admitted by ([1]) to integrating some 
linear system of PDEs (for further details, see [8]-|10j). 

However, with all importance and power of the traditional Lie approach, 
it does not provide all the answers to the mounting challenges of the mod- 
ern nonlinear physics. By this very reason there were numerous attempts 
of generalize Lie symmetries so that the generalized symmetries retain the 
most important features of Lie symmetries and allow for broader scope of 
applicability. The natural move in this direction would be to allow for the 
coefficients of infinitesimal generators of Lie symmetries to contain not only 
independent and dependent variables and their derivatives but integrals of 
dependent variables, as well. In this way, the so called nonlocal symmetries 
were introduced into mathematical physics. 

The concept of nonlocal symmetry of linear PDEs is well understood by 
now (see, e.g., [H]). However, this is not the case for nonlinear equations. 
The problem of developing regular methods for constructing nonlocal sym- 
metries of nonlinear PDEs is still waiting for its Sophus Lie. On the other 
hand, there is a number of results on nonlocal symmetries for specific equa- 
tions. One of the possible approaches to construction of nonlocal symmetries 
has been suggested by Bluman [12]- [H]. He put forward the concept of po- 
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tential symmetry, which is a special case of nonlocal symmetry. The basic 
idea of the method for constructing potential symmetries of PDEs can be 
formulated in the following way. Consider evolution equation 

Ut = f{t,X,U,Ui, . . . ,Un). (2) 

Suppose that it can be rewritten in the form of a conservation law 

^ (G{t, X, m) j = ^ (^F{t, X, U,Ui,..., Un-l)^ . (3) 

By force of we can introduce new dependent variable v = v{t,x) and 
rewrite equation ([1]) as follows 

Vx = G(t,X,u), Vt = F(t,X,U,Ui, . . . ,Un-l). (4) 

Now if system of two equations (jl]) admits Lie symmetry such that at least one 
of the coefficients of its infinitesimal operators depends on v = d~^G{t, x, u), 
then this symmetry is the nonlocal symmetry for the initial evolution equa- 
tion (I2l). Here is the inverse of dx, i.e., dx = dx = 1. This nonlocal 
symmetry is also called potential symmetry of ([2]). 

Pucci and Saccomandi proved that potential symmetries can be derived 
using non-classical symmetries of the ([2]). Recently, we established much 
stronger assertion by associating potential symmetries with classical (con- 
tact) symmetries [T7|. More precisely, we proved that any potential sym- 
metry of evolution equation ([2]) can be reduced to contact symmetry by a 
suitable nonlocal transformation of dependent and independent variables. As 
a consequence, one can obtain exhaustive description of potential symmetries 
of ([2]) through classification of contact symmetries of PDEs of the form ([2]). 

Some applications of potential symmetries to specific subclasses of ([2]) 
can be found in [TS]-[2B]. 

In the present paper we generalize results of [17J for system of evolu- 
tion equations ([1]) and prove that any potential symmetry of the system in 
question reduces to classical Lie symmetry under a suitable nonlocal trans- 
formation of dependent and independent variables (Sections 1,2). Next, we 
suggest in Section 3 a more general approach to constructing nonlocal sym- 
metries that goes far beyond of the concept of potential symmetries. This 
approach enables generating systems of evolution equations associated with 
a given system of the system ([1]), provided the latter admits non-trivial Lie 
symmetry. We give several examples of application of the approach in Section 
3. 
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2 Conservation law representation and clas- 
sical symmetries 

Definition 1 We say that system (J\) admits complete conservation law rep- 
resentation (CLR) if it can be written in the form 

^ (G{t, X, u) j = ^ (F{t, X, u, ui, . . . , u„_i)) . (5) 
Here u, F, G are m-component vectors. 

Definition 2 We say that system (J\) admits partial CLR if it can be written 
in the form 

^ (F{t, X, u, w)) = ^ (^G(t, X, u, ui, . . . , u„_i, w, wi, . . . , w„_i)) , 
Wf = H(t, X, u, ui, . . . , u„, w, wi, . . . , w„). (6) 

Here u, F, G and w, H are r-component and m — r-component vectors, re- 
spectively. 

Below we present theorems that provide exhaustive characterization of 
conservation law representability in terms of classical Lie symmetries. We 
give the detailed proof of the assertion regarding complete CLR, the case of 
partial CLR is handled in a similar way. 

Theorem 1 System (QP admits complete CLR if and only if it is invariant 
under m- dimensional commutative Lie algebra Cm = (ci, . . . , Cm), where 

m 

Ci = ^i{t, X, u)d^ + ^vl it, a;, u)9„, , (7) 
i=i 

and besides, 



rank 



m. 



(8) 
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Proof. Suppose system ([T]) admits CLR Introducing new m-component 
function 

V, = G(t,a;,u) (9) 
and eliminating u from ([5]) we get 

Vxt = -^(i{t,X,Vi,...,\n)y (10) 

Integrating the obtained system of PDEs with respect to x yields 

vt = f(t,a;, vi, . . . , v„). (11) 

Note that the integration constant w(t) is absorbed into the function v. 
Evidently, system (fTTj) is invariant under the commutative m-dimensional 
Lie algebra Cm = {dy'^, ■ ■ ■ , d^m). What is more, the coefficients of the basis 
elements of the algebra satisfy condition (jSj). 

Let us prove now the that the inverse assertion is also true. Suppose that 
([1]) admits Lie algebra Cm = (ei, . . . , e^), whose basis elements have the 
form ([7]) and satisfy ([8]). Then there is a change of variables (see, e.g. [8]) 

t = t, x = X(t, x, u), u = U(t, a;,u) 

reducing basis elements of Cm to the form Cj = 9^^*, i = 1, . . . ,m. In what 
follows we drop the bars. 

Now ([1]) necessarily takes the form 

= f(t,a;, ui,...,u„). (12) 

Differentiating (fT2|) with respect to x and making the (nonlocal) change of 
dependent variables = u, we finally get 

Vi = ^f(t,x, V, Vi, . . . , v„_i), 
which completes the proof. □ 

Note 1. The fact that symmetry operators, ei, . . ., are of specific form (jTl) 
is crucial for the whole procedure of reducing a system of evolution equations 
to a 'conserved' form ([5]). If a symmetry group generated by some operator 
does not preserve the temporal variable, t (which means that the coefficient 
of dt in Cj is non-zero for some i), then this operator cannot be reduced to 
the canonical form d„i and the reduction routine does not work. 
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Theorem 2 System (QP admits partial CLR if and only if it is invariant 
under r- dimensional commutative Lie algebra Cr = (ci, . . . , e^), where 



and besides, 



Ci = ^i{t, X, u, w)9^ + '7i a;, u, w)d^ 



^ CKt,x,u,w)d^j i = l, 



/a ^! ••• Ci' 



rank 



' Ir Sr 



cr \ 



(13) 



3 Potential symmetries 

Potential symmetries of system of evolution equations ([1]) appear in the same 
way as they do for a single evolution equation. For simplicity, we consider the 
case of complete CLR. By force of (jSD we can introduce the new dependent 
variable v, so that 



Vi = F{t, X, u, ui, . . . , u„_i), v^. = G(t, X, u). 

Note that v is nonlocal variable since v = d~^G{t, x, u). 
Suppose now that system (HM admits Lie symmetry 

t' = T{t, X, u, V, 6), X = X{t, X, u, V, 6), 
u' = U(t, X, u, V, 6), v' = V(t, X, u, V, 6), 



such that one of the derivatives 

dT dT dV d\ 



(14) 



(15) 



, m 



does not vanish identically. Rewriting group (|T5l) in terms of variables t, x, u 
and taking into account that v = d^^u yield the nonlocal symmetry of the 
initial system of evolution equations ([T]). This means, in particular, that sym- 
metry in question cannot be obtained within the Lie infinitesimal approach. 
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What we are going to prove, is that this symmetry can be derived by regular 
Lie approach if the later is combined with the nonlocal transformation of the 
dependent variables. 

Indeed, let system ([1]) admit complete CLR ([5]). In addition, we sup- 
pose that ([1]) possesses potential symmetry. Making the nonlocal change of 
dependant variables, u — > v, 

Va; = G(t,x, u), U = G(t,X, V^.), G(t,X,G(t,X,Y^)) (16) 

we rewrite dSD in the form flTUl) . As initial system ([T]) admits a potential 
symmetry, system (IT^ is invariant under the Lie transformation group of 
the form (ITSl) . 

Integrating (fTOl) with respect to x yields system of evolution equations 

vt = f(t,2;, vi, . . . , v„). (17) 

Next we rewrite Lie symmetry (1151) by eliminating u according to (I16p which 
yields 

t' = T(t,x,G(t,x,v^),\,9^, x' = x(t,x,G(t,x,\^),Y,9^, 

v' = v(t,a;,G(t,x,v,),v,^). (18) 

By construction. Lie transformation group (ITSl) maps the set of solutions of 
(fT7|) into itself. Consequently, (fTSi) is the Lie group of contact symmetries of 
system of evolution equations (fT7|) . 

It is a common knowledge that any contact symmetry of a system of PDEs 
boils down to the first prolongation of a classical symmetry [21]. Conse- 
quently, the derivatives of T, X, V with respect to the third argument vanish 
identically and we get 

t' = T{t,x,v,e), x' = X{t,x,v,9), V = V(t,x,v,^). (19) 

This group is nothing else than the standard Lie symmetry group of system 

The same assertion holds true for the case of partial CLR. 

Theorem 3 Let system of evolution equations (QP admit complete or partial 
CLR and be invariant under a potential symmetry. Then there exist a ( non- 
local) change of variables mapping ^ into another system of the form (Q]) 
so that the potential symmetry of ^ becomes the standard Lie symmetry of 
the transformed system. 



7 



This assertion is in fact a no-go theorem for potential symmetries of sys- 
tems of evolution equations. It states that the concept of potential symmetry 
does not produce essentially new symmetries. The system admitting poten- 
tial symmetry is equivalent to the one admitting standard Lie symmetry, 
which is the image of the potential symmetry in question. 

However, there is more to it. Theorem 3 imply the regular algorithm 
for group classification system of nonlinear evolution equations admitting 
nonlocal symmetries. Again, for the sake of simplicity, we consider the case 
of complete CLR. 

Indeed, let system of evolution equations ([1]) be invariant under (m + 
l)-dimensional Lie algebra Cm+i = (ci, Cm+i)- Here Ci, . . ., are 
commuting operators of the form ([7]) and their coefficients satisfy constraint 
([8]). Basis operator Cm+i is of the generic form 

m 

em+i = T{t, X, u)dt + ^i{t, X, u)d:, + ^ r]l {t, X, n)duj . 

i=i 

Making an appropriate change of variables we can reduce the operators 
Ci, . . ., Cm to the canonical forms, namely, Cj = d^t, i = l,...,m. Then 
system ([1]) necessarily takes the form f|T7j) . 

Let f|T9|) be Lie transformation group generated by the symmetry oper- 
ator Cm+i- Calculating the first prolongation of formulas f|T9l) we get the 
transformation rule for the first derivatives of v 

v; = W(t,a;,v,v,,0). (20) 

Now we differentiate (ITUl) with respect to x and make the following change 
of dependent variables equations, 

w = V,, (21) 

which yield 

Wi = — (f(t,x,w,...,w„_i)). (22) 

Formulas f|T9|) , fl20|) provide the image of the transformation group ([T9l) under 
the mapping ( l20l) . so that 

t' = T{t,x,v,e), x' = x{t,x,v,e), w'^ = w{t,x,v,w,e). (23) 

Here v = d~^w. 
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Consequently, if one of the derivatives, dT/dv\dX /dv\d'W/dv\ does 
not vanish identically, then (1231) is the nonlocal symmetry group of system 
of evolution equations fl22l) . 

The same line of reasonings applies to the case when system ([1]) admits 
partial CLR. 

We summarize the above speculations in the form of the multi-step algo- 
rithm for group classification of nonlocal symmetries of systems of evolutions 
equations associated with a given system of the form ([1]). 

Let system of evolution equations ([1]) be invariant under A^-dimensional 
Lie symmetry algebra Cn- For simplicity we formulate the algorithm for the 
case of complete CLR. 

Algorithm 1. Classification of potential symmetries of (QP 

1. Calculate inequivalent subalgebras Ai of the algebra C^. 

2. Select those subalgebras Ai, which contain commutative subalgebras 

of operators of the form ([7]). 

3. For each commutative subalgebra Aim perform change of variables re- 
ducing its basis elements to the canonical forms d^i, . . ., d^m and trans- 
form correspondingly initial system ([T]). 

4. Perform nonlocal transformation fl?Il) . 

5. Eliminate 'old' dependent variables v from f|T9l) in order to derive sym- 
metry group fl23l) of the transformed system of evolution equations fl22|) . 

6. Verify that there is, at least, one derivative from the list dT/dv"^, 
dX/dv^, dW/dv^ that does not vanish identically. If this is the case, 
then f l25]) is the nonlocal (potential) symmetry of fl22]) . 

The steps needed to implement the above algorithm for the case of system 
of evolution equations admitting partial CLR are the same, the only differ- 
ence is that intermediate formulas ffT^ - fl25]) are more cumbersome, since we 
need to distinguish between two sets of dependent variables u and w (see, 

©)• 

Note that by force of Theorems 1,2 any potential symmetry of equations 
of the form ([1]) can be obtained in the above described manner. 
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4 Some generalizations 



Denote the class of partial differential equations of the form ([T]) as Then 
any system of the form 

uj = /(t,x,ui, . . .u„), (24) 

(i) belongs to €„, and, (u) its image under nonlocal transformation u = v^; 
also belongs to Existence of such nonlocal transformation is in the core 
of our approach to classifying nonlocal symmetries of systems of evolution 
equations. 

It is not but natural to ask whether there are other types of nonlocal 
transformations of the class that can be utilized to generate nonlocal sym- 
metries. Remarkably, such nonlocal transformations do exist. Sokolov [25] 
put forward the idea of group approach to generating such transformations 
for a single evolution equation. It is straightforward to modify his approach 
to handle systems of evolution equations, as well. As an illustration, we 
consider system (12^ . It is invariant under the m-dimensional Lie algebra 
'^m = • • • , du^). The simplest set of (m + 2) functionally-independent 
invariants of the algebra Cm can be chosen as follows t, x, u]., . . ., m^. Now 
we define the transformation 

t T(t, X, U, Ux, Uj-xy • • ^(^5 -^y ^5 ^xxy • • •)) 

u = U(t,x,u,u^,u^.^., . . .) (25) 

so that T, X, U are invariants of the symmetry group of the system under 
study. In the case under consideration, we have T = t, X = x, \J = u^. As 
we established in Section 1, applying this transformation to any equation of 
the form yields system of evolution equations that belongs to (£„. What 
is more. Lie symmetry group of (12^ is mapped into symmetry group of the 
transformed system and some of the basis operators of the latter become 
nonlocal ones. 

Consider, as the next example system of evolution equations 

Uf = /()f:,x,U2,...u„), n > 3. (26) 

This system is invariant under the 2m-dimensional Lie algebra = {du^, 
. . ., du^, xdui, ■ . ., xdum). The simplest set of m + 2 functionally independent 
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first integrals reads as t,x, u].^, . . ., -u^. Consequently, change of variables 
f l25p takes the form 

t = t, X = X, V = u^^. (27) 

Note that we dropped the bars and replaced u with v. 
Transforming (126!) according to (1271) we get 

((^x^)^v)^ = f{t,X,V,Vi,. . .Vn-2) 

or, equivalently, 

(d-'f (vt - d'Jit, X, V, Vi, . . . V„_2)) = 0. 
Integrating twice yields 

Vt = dlf{t,X,V,Yi,...Yn-2)- (28) 

Note that integration constants w^(t)a; + w^(t) are absorbed by the function 

V. 

So that nonlocal transformation fl271) maps a subset of equations from 
into €n- Consequently, it can be used to generate nonlocal symmetries of the 
initial system fl2B]) . 

Let system (126!) be invariant under the Lie transformation group 

t' = T{t,x,u,e), X = x{t,x,u,e), u = v{t,x,u,e). (29) 

Computing the second prolongation of the above formulas we get the trans- 
formation law for the functions v = u^x, 

V = \(t,x,u,Ux,v,6). (30) 

Combining (!29l) and (130|) yields the symmetry group of system of evolution 
equations (!28l) . 

t' = T(t,x,u,9), X = X{t,x,u,9), v' = V(t, X, u, Ua;, V, 6*) (31) 
where u = ((9^^)^v are nonlocal variables. Now, if one of the derivatives 

dT dX dY dV 
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does not vanish identically, then fl3ip is the nonlocal symmetry group of 
system of evolution equations (125]) . 

It is important to emphasize that the symmetry algebra Cm is not obliged 
to be commuting. The necessary condition is that the corresponding trans- 
formation group has to preserve the temporal variable, t, i.e., basis elements 
of Lra have to be of the form 

m 

Q = ^{t,x,u)d^ + ^ r]^{t,x,u)d^j. (32) 

As an illustration, consider the following system of second-order evolution 
equations: 

. . / 7/1 7,''" \ 

«I = <f t,a;,^,...,^ , ^ = l,...,m. (33) 

This system is invariant under the 2m-dimensional Lie algebra = (<9m1, 
. . ., du^, u^du^, . . ., u^du^)- Note that the algebra is not commutative. 
The set of m + 2 invariants of the algebra £2™ can be chosen as follows, 

J. XX "'XX 

making the change of variables 

t = t, X = X, v' = ^,...,V^ = ^ 

we rewrite ( !33|) in the form 

-(a-^exp(a-V)) =exp(9-V)r(^,^,^\---,^'"), ^ = l,...,m. (34) 

Taking into account that the operators ^ and (9~^ commute, differentiating 
f lM|) with respect to x, and replacing v with we finally get 

w;j = <f(t,x,ti;^.,...,0 + ^r(t,x,w;^,...,0, 2=l,...,m. (35) 

The above system is obtained from the initial one through the change of 
dependent variables = 5"^ exp(w*), i = 1, . . . , m. Consequently, if system 



12 



(155]) admits symmetry fl29p . then system fl35p admits the following transfor- 
mation group: 

t' = T{t,x,u,6), X = X{t,x,u,6), w = W(t, X, u, w, ^) (36) 
with = exp(w*), i = 1, . . . ,m. Again, if one of the derivatives, 

dT dX dV dW 

du^' dv^' dv}' 

does not vanish identically, then ( 136|) is the nonlocal invariance group of 
system of evolution equations ( IMl) . 

The algorithm for group classification of nonlocal symmetries of system 
([Tj) suggested in the previous section yields those nonlocal symmetries which 
are potential, since the nonlocal transformation was chosen a priori. Allowing 
for a nonlocal transformation to be determined by symmetry group of the 
system under study, yields a more general algorithm for constructing nonlocal 
symmetries. 

Let system of evolution equations ([1]) be invariant under A^-dimensional 
Lie symmetry algebra C^- Then the following multi-step algorithm can used 
to construct nonlocal symmetries of ([T]). 

Algorithm 2. Classification of nonlocal symmetries of (J\) 

1. Calculate inequivalent subalgebras A4 of the algebra C^. 

2. Select those subalgebras A4, which contain basis elements ei, . . ., of 
the form fl32|l . 

3. For each Ai construct r + 2 functionally independent-invariants uj^{t, x, 

U, U^,...), Uj'^{t,X,U,U^,...), UJ^{t,X,U,U^,. . .), Uj''{t,X,U,U^,. ..) 

and make change of variables 

i=uj^, x = u^, 'ii = uj\ i = l,...,r. (37) 

4. Eliminate 'old' dependent variables u from fl37p in order to derive sym- 
metry group Q of the transformed system of evolution equations. 

5. Verify that there is, at least, one function from the list {cj*, u^, u^, . . ., 
cu^} that depends on for some 1 < i < r. If this is the case, then Q 
is the nonlocal symmetry of fl22l) . 
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5 Conclusion. 



One of the principal results of the paper is Theorem 3 from Section 2 stating 
that any potential symmetry of system of evolution equations ([1]) reduces 
to a Lie symmetry by an appropriate nonlocal transformation of dependent 
and independent variables. The nonlocal transformation in question is a 
superposition of the local change of variables 

i=t, x = X(t,x,u), u = U(t, x,u) (38) 

and of the nonlocal change of dependent variables 

V = u,. (39) 

The explicit form of transformations 0381) is defined by the Lie symmetry 
admitted by the corresponding system ([1]). 

we obtain as an by-product exhaustive characterization of systems ([1]), 
that can be represented in the form of conservation law(s), in terms of Lie 
symmetries preserving the temporal variable, t, 

t' = t, x' = X{t,x,u,e), u = U(t, x, u, ^) 

(see, Theorems 1,2). 

In Section 3, we generalize the above reasonings in order to obtain nonlo- 
cal symmetries which are not potential. The basic idea is replacing 0391) with 
a more general nonlocal transformation. This transformations is determined 
by invariants of Lie symmetry algebra of the system under study. 

We intend to devote one of our future publications to systematic study of 
nonlocal symmetries of systems of nonlinear evolution equations ([T]) within 
the framework of the approach developed in Section 3. 
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